Optimizing the discrete time quantum walk using a SU(2) coin 
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We present a generalized version of the discrete time quantum walk, using the SU (2) operation 
as the quantum coin. By varying the coin parameters, the quantum walk can be optimized for 
maximum variance subject to the functional form a 2 w iV 2 and the probability distribution in 
the position space can be biased. We also discuss the variation in measurement entropy with the 
variation of the parameters in the SU(2) coin. Exploiting this we show how quantum walk can be 
optimized for improving mixing time in an n-cycle and for quantum walk search. 



I. INTRODUCTION 

The discrete time quantum walk has a very similar 
structure to that of the classical random walk - a coin flip 
and a subsequent shift - but the behaviour is strikingly 
different because of quantum interference. The variance 
a 2 of the quantum walk is known to grow quadratically 
with the number of steps N, a 2 oc TV 2 , compared to the 
linear growth, a 2 oc N, for the classical random walk 
U S SHI • This has motivated the exploration for a new 
and improved quantum search algorithms, which under 
certain conditions are exponentially fast compared to the 
classical analog Q. Environmental effects on the quan- 
tum walk Q and the role of the quantum walk to speed 
up the physical process, such as the quantum phase tran- 
sition have been explored Q ■ Experimental implementa- 
tion of the quantum walk has been reported [8| and var- 
ious other schemes for a physical realization have been 
proposed @. 

The quantum walk of a particle initially in a 
symmetric superposition state \^i n ) using a single- 
variable parameter 9 in the unitary operator, Ug = 
cos(0) sin(0) 
sin(0) - cos(0) 
metric probability distribution in the position space. The 
change in the parameter 9 is known to affect the varia- 
tion in the variance, a 2 Q. It has been reported that 
obtaining a symmetric distribution depends largely on 
the initial state of the particle 0, EL G3 ■ 

In this paper, the discrete time quantum walk has been 
generalized using the SU(2) operator with three Caley- 
Klein parameters £, 9 and C as the quantum coin. We 
show that the variance can be varied by changing the pa- 
rameter 9, a 2 sa (1 — sin(0))iV 2 and the parameters £ and 
£ introduce asymmetry in the position-space probability 
distribution even if the initial state of the particle is in 
symmetric superposition. This asymmetry in the proba- 
bility distribution is similar to the distribution obtained 
for a walk on a particle initially in a non-symmetric su- 
perposition state. We discuss the variation of measure- 
ment entropy in position space with the three parame- 
ters. Thus, we also show that the quantum walk can be 
optimized for the maximum variance, for applications in 



, as quantum coin returns the sym- 



search algorithm, improving mixing time in an n-cycle 
or general graph and other processes using a generalized 
SU(2) quantum coin. The combination of the measure- 
ment entropy and three parameters in the SU(2) coin can 
be optimized to fit the physical system and for the rele- 
vant applications of the quantum walk on general graphs. 
This paper discuss the effect of SU(2) coin on quantum 
walk with particle initially in symmetric superposition 
state. The SU(2) coin will have a similar influence on 
a particle starting with other initial states but with an 
additional decrease in the variance by a small amount. 

The paper is organized as follows. SectionHTlintroduces 
to the discrete time quantum (Hadamard) walk. Sec- 
tion [Hi] discusses the generalized version of the quantum 
walk using the arbitrary three-parameter SU (2) quan- 
tum coin. The effect of three parameters on the variance 
of the quantum walker is discussed, and the functional 
dependence of the variance due to parameter 9 is shown. 
The variation of the entropy of the measurement in po- 
sition space after implementing the quantum walk using 
different values of 9 is discussed in Sec. |IV] Section IVl 
and I VII discuss optimization of the mixing time of the 
quantum walker on the n— cycle and the search using a 
quantum walk. Section I VIII concludes with a summary. 



II. HADAMARD WALK 



To define the one-dimensional discrete time quantum 
(Hadamard) walk we require the coin Hilbert space Ti. c 
and the position Hilbert space Ti. p . The TL C is spanned 
by the internal (basis) state of the particle, |0) and |1), 
and the TL P is spanned by the basis state \ipi), i G Z. 
The total system is then in the space TC — 7i c eg) Tt p . 
To implement the simplest version of the quantum walk, 
known as the Hadamard walk, the particle at origin in 
one of the basis state is evolved into the superposition of 
the Ti. c with equal probability, by applying the Hadamard 
1 1 



operation, H = 
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such that, 
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(#®1)(|0>®|V>o>) = ^[|0) + |1>]®K>o> 
(ff <8> 1)(|1) ® |Vo)) = ^=[|0) - |1)] (8 |Vo>- (1) 

The is then followed by the conditional shift operation 
5: conditioned on the internal state being |0) (|1)) the 
particle moves to the left (right), 

S = |0><0| ® £|^_i><V»i| + |1)<1| ® (2) 

The operation S evolves the particle into the superposi- 
tion in position space. Therefore, each step of quantum 
(Hadamard) walk is composed of an application of H and 
a subsequent S operator to spatially entangle TL C and TL p . 
The process of W = S(H <g> 1) is iterated without resort- 
ing to the intermediate measurements to realize a large 
number of steps of the quantum walk. After the first 
two steps of implementation of W, the probability dis- 
tribution starts to differ from the classical distribution. 
The probability amplitude distribution arising from the 
iterated application of W is significantly different from 
the distribution of the classical walk. The particle with 
initial coin state |0) (|1)) drifts to the right (left). This 
asymmetry arises from the fact that the Hadamard oper- 
ation treats the two states |0) and |1) differently, multi- 
plies the phase by —1 only in case of state |0). To obtain 
left-right symmetry in the probability distribution, (b) 
in Fig. (jT| b), one needs to start the walk with the par- 
tilce in the symmetric superposition state of the coin, 



|*<n> 



V2 



[|0)+i|l)]®hfo>. 



III. GENERALIZED DISCRETE TIME 
QUANTUM WALK 

The coin toss operation in general can be written as an 
arbitrary three parameter SU (2) operator of the form, 



cos(0) 
sin(0) 



3<S 



sin(0) 
cos(0) 



(3) 



the Hadamard operator, H = J7o.f ,o- By replacing the 
Hadamard coin with an operator Ugg^, we obtain the 
generalized quantum walk. For the analysis of the gener- 
alized quantum walk we consider the symmetric superpo- 
sition state of the particle at the origin. By varying the 
parameter £ and £ the results obtained for walker starting 
with one of the basis (or other nonsymmetric superposi- 
tion) state can be reproduced. A particle at origin in a 
symmetric superposition state |\I/m), when subjected to 
a subsequent iteration of W^ t e,c = X ® ■"■) imple- 

ments a generalized discrete time quantum walk on a line. 
Consider an implementation of W^g t t, which evolves the 
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(b) = (0°, 45°, 0°) 

(c) = (0°, 60°, 75°) 
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FIG. 1: The spread of probability distribution for different 
value of 9 using operator Uo t e,o, is wider for (a) = (0, j^, 0) 
than compared to (b) = (0, ?,0). Biasing the walk using ( 
shifts the distribution to right, (c) = (0, -|, and £ shifts to 
the left, (d) = (ff , f ,0). The distribution is for 100 steps. 



walker to, 

Wt AC \* in ) = 4=[(e^ cos(0) +ie« sin(0)) |0>|V>-i> 

+ (e~ lC sin(0) - ie-^ cos(0)) (4) 

If £ — C, Eq. (01 has left-right symmetry in the position 
probability distribution, but not otherwise. We thus find 
that the generalized SU{2) operator as a quantum coin 
can bias a quantum walker in spite of the symmetry of 
initial state of the particle. We return to this point below. 

It is instructive to consider the extreme values of the 
parameters in the U^ t g^. If £ = = f = 0, ?7o,o,o = Z, 
the Pauli Z operation, then W^ex = & an( i * ne two 
superposition states, |0) and |1), move away from each 
other without any diffusion and interference having high 



N 2 . On the other hand, if 



then U ,z.,o=X, 



the Pauli X operation, then the two states cross each 
other going back and forth, thereby remaining close to 
position i = and hence giving very low cr 2 w 0. These 
two extreme case are not of much importance, but they 
define the limits of the behavior. Intermediate values of 
the between these extremes show intermediate drifts 
and quantum interference. In Fig. |T]) we show the sym- 
metric distribution of quantum walk at different values 
of by numerically evolving the density matrix. Fig. @ 
shows the variation of a with increase in for quantum 
walk of different number steps with the operator, J7o,e,o- 
The change in the variance for different value of the is 
attributed to the change in the value of Cg , a constant for 
a given 0, a 2 = CeN 2 , Fig. Therefore, starting from 
the Hadamard walk (0 = j; £ = ( = 0), the variance can 
be increased (0 < f ) or decreased (0 > f ) respectively. 
In the analysis of Hadamard walk on the line in Q , it 
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FIG. 2: A comparison of variation of a with 8 for different 
number of steps of walk using operator Uo,e,o using numerical 
integration. 



FIG. 3: The probability distribution obtained using using Eq. 
(f5]) for different value of 9. The distribution is for 100 steps. 



is shown that after N steps, the probability distributed is 
spread over the interval ^] and shrink quickly out- 
side this region. The moments have been calculated for 
asymptotically large number of steps N and the variance 

is shown to vary as a 2 (N) = (l — -1=J N 2 [3|. 

The expression for the variance of the quantum walk 
using £/o,0,o as a quantum coin can be derived by us- 
ing the approximate analytical function for the proba- 
bility distribution P(i) that fit the envelop of the quan- 
tum walk distribution obtained from the numerical inte- 
gration technique for different values of 9. For a quan- 
tum walk using t/0,0,0 as quantum coin, after N steps 
the probability distribution is spread over the interval 
(-iVcos(0), Ncos(0)) [I] . This is also verified by the an- 
alyzing the distribution obtained using the numerical in- 
tegration technique. By assuming the value of the prob- 
ability to be beyond |iVcos(0)|, the function that fits 
the probability distribution envelop is, 



ATcos(0) 



Ncos(B) 



[1+ COS 2 (20)] 



K{6)( 



N 



P(i)di ■ 



-di w 1, (5) 



where, K(B) = cos(0)[l + cos 2 (2(9)] [1 + sin(0)] [H. 
Fig. shows the probability distribution obtained by 
using the Eq. (©. The interval (-jVcos(0), iVcos(0)) 
can be parametrised as a function of <j), i = f(<f>) — 
N cos(0) sin(0) where <f> range from — -| to J. For a walk 
with coin C/o,e,Oj the mean of the distribution is zero and 
hence the variance can be analytically obtained by eval- 
uating, 



/N cos(S) 
P(i)i 2 di 
-jvcos(e) 



P(/W)(/W>))VW#-(6) 




FIG. 4: Variation of Cg when rj = |£ — f | = 0° from numerical 
integration to the function (1 — sin(0)) to which it fits. The 
effect of maximum biasing, r\ — 90° on Cg is also shown and 
its effect is very small. 



(l + COS 2 (20)) eg(e)(sln2w _ l) ^^^y 

Vn 



x (7Vcos(0) cos(0)) d<f) = N 2 (l - sin(0)).(7) 



cr 2 = C N 2 



(l-sin(0))iV 2 . 



(8) 



We also verify from the results obtained through numer- 
ical integration that Cg = (1 — sin(0)), Fig. 

Setting £ 7^ £ in U^g^ introduces asymmetry, biasing 
the walker. Positive £ contributes for constructive inter- 
ference towards right and destructive interference to the 
left, whereas vice versa for £. The inverse effect can be 
noticed when the £ and £ are negative. As noted above, 



4 



for £ = C, the evolution will again lead to the symmetric 
probability distribution. Apart from a global phase, one 
can show that the coin operator 



o,e,C-€- 



(9) 



In Fig. ([T|) we show the biasing effect for (£, 6*, £) = 
(0°, 60°, 75°) and for (75°, 60°, 0°). The biasing does not 
alter the width of the distribution in the position space 
but probability goes down as a function of cos(?7) on one 
side and up as a function of sin(?7) on the other side. 
Where r] = |£ — CI ■ The mean value i of the distribution, 
which is zero for f7o,e,Oi attains some finite value with 
non-vanishing rj, this contributes for an additional term 
inEq. ©, 



r jVcos(e) 
<j- - / P{i){i-ifdi. 

-Ncos(8) 



(10) 



this contributes to a small decrease in the variance of the 
biased quantum walker, Fig. 

It is understood that, obtaining symmetric distribu- 
tion depends largely on the initial state of the particle 
and this has also been discussed in 0, 0, El E| ■ But 
using f/{,0,£ as coin operator, and examining the walk 
evolution shows how non-vanishing £ and C introduce 
bias. For example, the position probability distribution 
in Eq. ((4J corresponding to the left and right positions 
are j[l± sin(20) sin(£ — C)], which would be equal, and 
lead to a symmetric distribution, if and only if £ = £• 
The evolution of the state after n steps, [Wg 9,t] n | 1 ^*n) i s 
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FIG. 5: Variation of entropy of measurement H(i) with 9 for 
different number of steps. The decrease in H(i) is not drastic 
till 6 is close to or -|. 



position probability distribution pi obtained by tracing 
over the coin basis: 



H{i) = -^Pilogpi 



(14) 



|*(n)>= Y, OW|0)|V>„ 



\m m )) (ii) 



m——n 



and proceeds according to the iterative relations, 

A m , n = e li cosM m _i,„_i + e lC sin^S m _i, n _i(12a) 
B m ,n = e~ lC cos0A m _i >n _i - e~' iC sin^S m _i, n _i(12b) 

A little algebra reveals that the solutions A m , n and 
£? m ,„ to Eqs. (fT2"|) can be decoupled (after the initial 
step) and shown to satisfy 

A m , n +i - 4n,»-i = cos6»(e' c yl m _i^ - e ie A m+ i,„()13a) 
B m , n +i — -Bm.,n-i = cos 6{e^B m -i >n — B m+ \ tn \13b) 

For spatial symmetry from an initially symmetric super- 
position, the walk should be invariant under an exchange 
of labels <-> 1, and hence should evolve A m ^ n and B m ^ n 
alike (as in the Hadamard walk [Oj])- From Eq. (|13[) . we 
see that this happens if and only if £- = C- 



IV. ENTROPY OF MEASUREMENT 

As an alternative measure of position fluctuation to 
variance, we consider the Shannon entropy of the walker 



The quantum walk with a Hadamard coin toss, J7o,i,0) 
has the maximum uncertainty associated with the prob- 
ability distribution and hence the measurement entropy 
is maximum. For £ = C = and low 9, operator Uo,e,o 
is almost a Pauli Z operation, leading to localization of 
walker at ±JV. At 6 close to |, with £ = C = 0, U ap- 
proaches Pauli X operation, leading to localization close 
to the origin, and again, low entropy. However, as 9 ap- 
proaches ?, the splitting of amplitude in position space 
increases towards the maximum. The resulting enhanced 
diffusion is reflected in the relatively large entropy at 2L, 
as seen in Fig. §5§. Fig. © is the measurement entropy 
with variation of 9 in the coin Uq^ } o for different num- 
ber of steps of quantum walk. The decrease in entropy 
from the maximum by changing 9 on either side of j is 
not drastic untill the 9 is close to or Therefore for 
many practical purposes, the small entropy can be com- 
pensated for by the relatively large Cg , and hence a 2 . For 
many other purposes, such as mixing of quantum walk 
on an n-cycle Cayley graph, it is ideal to adopt a lower 
value of 9. The effect of £ and £ on the measurement 
entropy is of very small magnitude. These parameters 
do not affect the spread of the distribution and the vari- 
ation in the height reduces the entropy by a very small 
fraction. 
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FIG. 6: A comparison of mixing time M of the probability 
distribution of a quantum walker on a n-cycle for different 
value of 6 using coin operation I/o,e,o> where n, the number 
of position, is 101. Mixing is faster for lower value of 9. The 
distribution is for 200 cycles. 



maximum spread in the cycle during each cycle. Maxi- 
mum spread during each cycle distributes the probability 
over the cycle faster and this would optimize the mixing 
time. Thus optimizing mixing time with lower value of 
9 can in general be applied to most of the finite graphs. 
For optimal mixing time, it turns out to be ideal to fix 
£ — C m ^£,0,0 since biasing impairs a proper mixing. 
Fig. ^ is the time averaged probability distribution of 
a quantum walk on an n-cycle graph after n log n time 
where n is 101. It can be seen that the variation of the 
probability distribution over the position space is least 
for 9 = 15° compared to 6 = 45° and 9 = 75°. 



VI. QUANTUM WALK SEARCH 

A fast and wide spread defines the effect of the search 
algorithm. For the basic algorithm using discrete time 
quantum walk, two quantum coins are defined, one for 
a marked vertex and the other for an unmarked vertex. 
The three parameter of the SU (2) quantum coin can be 
exploited for an optimal search. 



V. QUANTUM WALK ON THE n- CYCLE AND 
MIXING TIME 

The n-cycle is the simplest finite Cayley graph with n 
vertices. This example has most of the features of the 
walks on the general graphs. The classical random walk 
approaches a stationary distribution independent of its 
initial state on a finite graph. Unitary (i.e., non- noisy) 
quantum walk, does not converge to any stationary dis- 
tribution. But by defining a time-averaged distribution, 

t=0 

obtained by uniformly picking a random time t between 
and (T— 1), and evolving for t time steps and measuring 
to see which vertex it is at, a convergence in the prob- 
ability distribution can be seen even in quantum case. 
It has been shown that the quantum walk on an n-cycle 
mixes in time M = O(nlogn), quadratically faster than 
the classical case which is 0(n) [HI]. From Eq. ((6]) we 
know that the quantum walk can be optimized for max- 
imum variance and wide spread in position space, be- 
tween (— Ncos(9), N cos(9)) after TV steps. For a walk 
on an n-cycle, choosing 9 slightly above would give the 



VII. CONCLUSION 

In this paper we have generalized the Hadamard walk 
to a general discrete time quantum walk with a SU (2) 
coin. We conclude that the variance of quantum walk can 
be optimized by choosing low 9 without loosing much on 
measurement entropy. The parameters £ and C intro- 
duce asymmetry in the position space probability distri- 
bution starting even from an initial symmetric superpo- 
sition state. This asymmetry in the probability distri- 
bution is similar to the distribution obtained for a walk 
on a particle initially in a non-symmetric superposition 
state. Optimization of quantum search and mixing time 
on an n-cycle using low 9 is possible. The combination of 
the parameters of the SU(2) coin and the measurement 
entropy can be optimized to fit the physical system and 
for the relevant applications of the quantum walk on a 
general graph. 
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